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Abstract-TYas paper deals with the problem of nonlinear dynamics 
of complex system with phase control combining phase-locked 
loop and automatic gain control loop. The behavior of the 
examined system is described by nonlinear two- and four- 
dimensional sets of differential equations with periodical 
nonlinearities. Stability of synchronous mode, bifurcations 
determining boundaries of domain with quality different behavior 
of the considered system, and oscillatory modes arising in domain 
where synchronous mode is unstable are studied. Results are 
presented using two-parameter bifurcation diagrams, phase 
portraits of attractors, time realizations of oscillations, and 
Poincare maps. 
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I. INTRODUCTION 

At present, the systems with phase control are widely used. 
As known [Shakhgil'dyan and Lyakhovkin, 1972; Lindsey, 
1972], such systems are traditionally intended to provide for 
and maintain the synchronous state, when the phase difference 
of reference and controlled signals becomes constant or, 
equivalently, the frequency difference of these signals is equal 
to zero. The systems may be also operating in non- 
synchronous modes with variable phase and frequency errors. 
The use of such modes opens wide possibilities for some 
nontraditional engineering and technological applications of 
phase control systems (generation of complex periodic and 
chaotic signals, data transmission and processing, etc. 
[Dmitriev and Shirokov, 2004; Dmitriev, Kletsov and Kuzmin, 
2009]). 

Among the studies of the systems with phase control, the 
study of dynamical modes, bifurcation, and transitions to 
chaotic behavior is extremely interesting. In this paper we 
investigate dynamical behavior of complex system with phase 
control combining phase-locked loop (PLL) and automatic 
gain control loop (AGCL). Different versions of such systems 
are of interest because they represent a circuit implementation 
of optimal algorithms for tracking estimation of variable 
parameters (phase angle 9(t) and amplitude A(t) of receiving 
signal [Tikhonov and Kulman, 1975; Kulman, Zheronkina, 
1969]). 

Combined PLL and AGCL operate on the basic of 
synchronization between receiving (estimated) and reference 
signals. For that, PLL performs coincidence of the signals 
frequencies, and AGCL regulates amplification factor in the 
phase control loop so as to reduce the influence of input signal 
amplitude alternation on PLL operation. This paper is mainly 
devoted to investigation of the characteristics of the system's 



dynamical behavior that are due to application of first- and 
second-order low-frequency filters (LFFs) in control circuits, 
and to coupling via phase and gain control circuits. 

II. THE SYSTEM MODELS UNDER CONSIDERATION 

Equations describing the dynamics of the considered 
combined PLL and AGCL can be derived from the equations 
for estimates 9* and A* of input signal parameters 9 and A 
obtained in [Tikhonov and Kulman, 1975; Kulman, 
Zheronkina, 1969]. These equations can be represented for 
mismatch (p =9(t)-9*(t) and amplitude ratio x=A*/A in the 
operator form (p=d/dt) as follows [Ponomarenko, 1986] 



Ptp/k, =y-K 1 (p)xsincp, j 
x = K 2 (p)[G (x) +k 2 (coscp+ ax sincp)] .1 



(1) 



In Equation (1), ki and k 2 are the amplification factors of 
control circuits; y is the initial frequency mismatch; Go(x) 
depends on the gain control circuit structure and is expressed 
in the forms of Go(x)=5o-k 2 x and Go(x)=o/x-k 2 x [Tikhonov 
and Kulman, 1975], 5o=A(/A (A is average input signal 
amplitude); a is the gain control circuit parameter; Ki(p) and 
K 2 (p) are the transfer functions of the LLFs in the PLL and 
AGCL control circuits; and a is the factor of coupling via 
control circuits. By complexity of Equations (1), we consider 
amplitude A as constant parameter having arbitrary positive 
values. 

The types of LFFs are determined by the models using for 
description of dynamics of parameters 9 and A. Consider the 
simplest first-order filters with transfer functions Ki(p)=l, 
K 2 (p)=l/(1+T p), where T is a time constant. Equation (1) 
may be written in this case as 



dcp/dx = y-x sincp, 

dx/dx = X(G(x) + cos(p + axsin(p), 



(2) 



where T=k\t is dimensionless time, G(x) is represented by 
functions: G(x)=5~Px and G(x)=a/x -|3x, p=l+l/k 2 , 5=5(/k 2 , 
X=k 2 /(Toki). System (2) is a dynamical system determined on 
cylindrical phase surface U ={(p(mod27i),x}. 

Now consider the second-order filter in the PLL and first- 
order filter in the AGCL with transfer functions 
K 1 (p)=l/[l+(T 1 +T 2 )p+T 1 T 2 p 2 ], K 2 (p)=l/(1+T p), where Tj and 
T 2 are the time constants. In this case, equation (1) are written 

as 
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d(p/dx = y, dy/dx = z, 
|idz/dx = y-xsincp-y-EjZ, 
s 2 dx/dx = G(x) + cos(p + axsin(p, 



(3) 



where s 1 =(T 1 +T 2 )k 1 , £ 2 =T k 1 /k2, |i=T 1 T 2 k 1 . System (3) is a 
dynamical system with four-dimensional cylindrical phase 
space U={cp(mod27i),y,z,x}. 

Further, we will consider the Models (2) and (3) for the 
values of variable x>0. Since Models (2) and (3) are nonlinear, 
its nonlocal study encounters serious difficulties. Therefore, 
we apply qualitative-numerical methods of nonlinear dynamics 
and computer simulation. 

III. DYNAMICAL STATES OF MODEL 

At first, we discuss the problem of local stability of 
synchronous mode. The equilibrium states of System (2) 
situated in the half phase cylinder x>0 are determined from the 
equations 



y - xsincp = 0, G(x) + coscp + axsinq) = 0. 



(4) 



Consider Equations (4) as function G(x)=5~Px. In this case, 
from Equations (4) and the form of G(x), it follows that, for 
0<a<P and the values of y,5,P,aeCo where 
Co={yi(5,P,a)<y<y 2 (5,P,a)}, system (2) has two equilibrium 
states Ai((pi,Xi) and A 2 ((p 2 ,x 2 ). The coordinates cpi, x l5 cp 2 , and 
x 2 are defined from Equation (4). For a>P, the domain 
Co={y>yi(8,P,a)}. If a=0, the boundaries of domain C 



Yi(5,P) = - 



(Vs 2 +8 + 35) yl\6^JW+8^Sf 



16p 



y 2 (5,P) = -y,(8,P). 



(5) 



Bifurcation values y=y, and y=y 2 corresponding to merging 
of equilibrium states Ai((pi,Xi), and A 2 (cp 2 ,x 2 ) are defined from 
equality 

yi(5,p,a)=sin(p m2 (5+cos(p m2 )/(P-asin(p m2 ), i=l,2. (6) 

In Equality (6) cpm2 is the solution of the equation 

Pcos(p(8+cos(p)-sin 3 (p(P-asin(p)=0. (7) 

The values of (pm2 are satisfied following conditions: 
-7t/2<(pml<0, 0«pm2<7i/2. Fig. 1 shows qualitative disposition 
of bifurcation curves y, and y 2 , and domain Co on the (y,8) 
plane for 0<a<P (Fig. la) and a>P (Fig. lb). 

5 |8 




a 




a 



y 



Fig. 1 Disposition of the domain C : a - 0<a<P, b - a>P 

Now consider Equations (4), if function G(x)=o/x-Px. 
From (4) and the form of G(x), it follows that, for 0<a<P and 
y,5,p,aeD , where D o ={a>max(a-(y,p,a),0, a + (y,p,a))}, 



System (2) has also two equilibrium states Al(cpi,Xi) and 
Al((p 2 ,x 2 ). For a>P, the domain is D ={cj>max(cj-(y,p,a),0)}. 
If parameter a=0, a + (y,P,a)=a-(-y,P,a). Bifurcation values 
g-g- and cj=cj + corresponding to merging of equilibrium 
states Ai(cpi,xi) and A 2 ((p 2 ,x 2 ) are defined by 



J+ _ y 2 (P-asin(p m )-ysin(p m cos(p+ 
sin2 <Pm 



y (P - a sin cp m ) - y sin y m cos cp„ 
sin 2 (p- 



o = - 



(8) 



In (8) 9m and ^ ( ^ e(0,7i/2), G (- n /2,0)) are the 
solutions of the equation 



y (2p-asin(p)coscp -sincp = 0. 



(9) 



Disposition of bifurcation curves cj=cj- and g=g+, and 
domain D on the (y,a) plane is qualitatively similar to 
disposition of curves y=y ; and y=y 2 on the (y,5) plane given in 
Fig. 1. 

The characteristic equation for system (2) is 

X 2 + qx + r = 0. (10) 
In (10), if function G(x) takes the form G(x)=5~Px, then 

q = A,(P-asin(Pj 2 ) + x i 2 cos( Pi 2 ' 
r = Mxj 2 cos(pj 2 (P-asincpj 2 )+ \ (11) 

+ (aXj 2 cos (pj 2 - sin (pj 2 ) - sin (pj 2 ) . 

If function G(x) takes the form G(x)=g/x~Px, then in 
Equation (10) 

q = X{gI x 2 2 + P~asin(pj 2 ) + Xj 2 coscpj 2 , 
r = X(Xj 2 coscpj 2 (<t/x 2 2 +P-asin(pj 2 )~ \ (12) 
-sincpj 2 (sin(pj 2 ~ ax i 2 cos( Pi 2))- 

By investigating the roots of the characteristic Equation 
(10), we find that the equilibrium state Ai((pi,Xi) are stable, 
whereas the equilibrium state A 2 (cp 2 ,x 2 ) is unstable of saddle 
type. 

The equilibrium state A, corresponds to synchronous mode 
of combined PPL and AGCL. If this mode is realized in this 
system, then the system can carry out tracking of estimated 
input signal's parameters. The parameter region where such 
mode exists is limited by the values of parameters y, 5, P, a, 
and g those belong to domains Co and D . Values cpj and Xi 
characterize the accuracy with which the parameters of input 
signal are estimated. 

Further, we consider nonlocal dynamic processes evolving 
in Model (2) when parameters y, 5, and g vary and remaining 
parameters are fixed. The results of qualitative-numerical 
investigation of Model (2) for G(x)=5~Px and for 
G(x)=o/x~Px reveal qualitatively similar dynamical states, 
bifurcation, parametrical portraits on the (y,5) plane and (y,o) 
plane, and corresponding phase portraits. Fig. 2 shows 
qualitative disposition of bifurcation curves on the (y,5) plane 
within the domain Co obtained for System (2), if function G(x) 
takes the form G(x)=5~Px and 0<a<p. 
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parameters belong to domain Bs, the combined PLL and 
AGCL is in synchronous mode irrespective of the initial 
conditions. 



Fig. 2 Bifurcation curve and domains 
of dynamical modes 

In Fig. 2 lines y~ and y* correspond to saddle A 2 

separatrix Si and S2 loop 110 of rotatory type that encompasses 
the phase cylinder U0. By investigating of saddle value: 
a s = -x 2 (y,5,p,a)cos(p 2 (y,5,p,a) -X$ + 

+ |iasin(p 2 (y,5,p,a), (13) 

we find that the loop 110, that arises at y= y~ , is stable 
(saddle value rjs<0), whereas the loop 110, that forms at y= y+ , 
is stable at the part of the curve y+ below the point M (where 
<js<0) and unstable at the part of the curve above the point 

M (where rjs>0). The point M corresponds to 5=5s(y, P, a, X), 
that defined from the equation rjs(y,5,P,a,X)=0. 

Therefore, upon passing throw the curve y~ , if y is 
decreased, a stable rotatory type (2n -periodic in cp) limit cycle 
L0 appears in the half phase cylinder x>0. When, as the result 
of y increasing, the System (2) crosses line y+ and 5<8s, a 
stable rotatory type limit cycle L0 also appears in the half 
phase cylinder x>0. Upon passing throw the curve y+ , as y is 

decreased and 5>5s, an unstable rotatory type limit cycle TO 
appears in the half phase cylinder x>0. 

Line yc starting from the point M corresponds to double 
rotatory type limit cycle in the half phase cylinder x>0. If, as 
the result of increasing y, line yc is crossed, a stable rotational 
type limit cycle L0 and an unstable rotatory type limit cycle TO 
are born in the half phase cylinder x>0. Upon passing through 
the line yc, as y is decreased, limit cycles L0 and TO merge and 
disappear. 

Limit cycle L0 is associated with the asynchronous mode 
of the combined PLL and AGCL such that phase error (p 
rotates and amplitude ratio x periodically oscillate about 
certain mean value. 

Bifurcation curves yi, y~ , ,yc, and y2 given at Fig. 2 

identify various parameter domains possessing qualitatively 
different dynamics of the Model (2). Fig. 3 shows the phase 
portraits of the System (2) for dynamic mode domains 
represented in Fig. 2. For the parameters from domain 
Bs={ y~ <y<min( ,yc)} there are no limit cycles and 
equilibrium state A t of System (2) is globally asymptotically 
stable in the half phase cylinder x>0. The phase trajectories of 
System (2) converge to A[ independently of initial state of the 
system (Fig. 3a). Therefore, when the values of the system 




Fig. 3 Phase portraits of model (2) 

For the parameters from domains B p ={yi<y< y^ } and 
={ y+<y<y2} a stable equilibrium state A] and stable limit 

cycle L0 are simultaneously exist on the phase cylinder U0 
(Fig. 3b). Depending on initial state, the trajectories of System 
(2) converge either to equilibrium state A 1 or to limit cycle L0. 
Therefore, depending on initial conditions, synchronous or 
asynchronous mode corresponding to these attractors on the 
phase cylinder U0 develops in the combined PLL and AGCL. 
The regions of attraction of equilibrium state Ai and limit 
cycles L0 are bounded by of saddle A 2 separatrix on the phase 
cylinder U0. For the parameters from domain Bc={yc<y<y^ } 

there are equilibrium state Aj and limit cycles L0 and TO on 
the phase cylinder U0 (Fig. 3c). The region II(A1) of attraction 
of equilibrium state Ai is bounded by unstable limit cycle TO. 
Therefore, if the initial conditions belong to domain located 
above limit cycle TO, then the synchronous mode is realized in 
combined PLL and AGCL. 

From the results obtained in framework of Model (2) we 
infer that, coupling via feedback signals in combined PLL and 
AGCL with the first-order filters make possible the appearance 
of asynchronous mode of rotational limit cycle, which is not 
possible in partial PLL and AGCL. 

IV. PECULIARITIES OF MODEL (3) DYNAMICAL BEHAVIOR 

If second-order LFF is used in the PLL, the dynamics of 
combined PLL and AGCL becomes drastically more complex. 
In addition to the synchronous mode and the periodic 
asynchronous mode, complex-periodic and chaotic 
asynchronous modes appear. Apart from this, the synchronous 
mode may lose stability and a transition to a periodic quasi- 
synchronous mode characterized by an oscillatory-type limit 
cycle may occur, the quasi-synchronous mode may become 
chaotic. Besides, the system may exhibits quasi-synchronous 
and asynchronous modes of quasi-periodic types 
corresponding to oscillatory and rotational two-dimension tori 
in the phase space. 

First, let us analyze the stability of the combined PLL and 
AGCL synchronous mode. System (3) with parameters 
y,5,P,aeCo and y,a,P,aeD has two equilibrium states 
Ai(<pi,0,0,xi) and A 2 ((p 2 ,0,0,x 2 ) located within the range x>0 of 
phase space U; coordinates (pi,Xi,cp 2 , and x 2 are determined 
from (4). Equilibrium state A[ may be either stable or unstable, 
while equilibrium state A 2 is unstable of saddle type. 

The conditions under which equilibrium state Aj is stable 
are determined by the roots of characteristic 

Equation: 
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X 4 + aa 3 + a 2 X 2 + a 3 X + a 4 = 0, 



(14) 



where 



a; = e, / |i - (G ' (x n ) + a sin cp n ) / 8 2 , 

a 2 = (l-8,(G'(x n ) + asin(p n )/s 2 )/ |i, 

a 3 =(x n cos(p n -(G'(x n ) + asin(p n )/£ 2 ))/|i, 

a 4 =-(sin 2 cp n +x n G'(x n )cos(p n )/(iae 2 ). 



(15) 



Applying the Routh-Hurwits criterion to the Equation (14), 
we obtained that equilibrium state A ; is stable for the values of 
parameters such that the following inequalities hold: 

ai,a 2 ,a 3 ,a 4 >0, aia 2 -a 3 >0, a 3 (aia2-a 3 )-ai 2 a4>0. (16) 

If Condition (16) are satisfied, the combined PLL and 
AGCL described by Model (3) have a synchronous mode 
corresponding to equilibrium state Ai((pi,0,0,xi). 

The domain of parameters CsO where Conditions (16) are 
satisfied corresponds to the region where the synchronous 
mode persists. Figs. 4 and 5 show boundary of the domains Co 
and D where the equilibrium states exist (curve y t ) and 
boundary of domain CsO where equilibrium state Ai is stable 
(curve ys), for G(x)=5~Px (Fig. 4) and G(x)=a/x-px (Fig. 5) 
correspondently. 
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Fig. 4 Bifurcation diagram (|j.,y) for G(x)=5-Px 




Fig. 5 Bifurcation diagram (n,y) for G(x)=cf/x-Px 

The curves yi and ys, and other bifurcation curves are 
plotted on plane (|i,y) for the values of parameters a=2, P=l.l, 
8i=l, £ 2 =2, 5=1.25, o=0.5. The domain CsO is located between 
the lines ys and yi. 

The domain CsO is divided by the curve yp into domains Cs: 
{max (yi,yp)<y<ys) and Cs 1 :{y 1 <y<min(yp, ys)}. The curve yp 
corresponds to the bifurcation of saddle-focus A 2 separatrix 



loop rotatory type (Rexi, 2 <0, lmxi, 2 ^0, %3<0, %4>0, where 
% 1,2,3,4 are the roots of the characteristic equation (14) for the 
equilibrium state A 2 ). For the parameters from domain C s , the 
equilibrium state Ai is only attractor in the phase space U. The 
trajectories in the phase space converge to Ai independently of 
initial state of the system. Therefore, the domain C s is a locking 
region of the combined PLL and AGCL. When the values of 
the Model (3) parameters belong to domain C s , the system is in 
the synchronous mode irrespective of the initial conditions. 

Now let us consider qualitative change of the system's 
behavior that occurs owing to variations in y and |i during the 
exit from domain C s . When, as a result of growing y or |i, the 
system crosses boundary ys, Conditions (16) are violated and 
System (3) exhibits the Andronov-Hopf supercritical 
bifurcation. The latter is related to a solution to characteristic 
Equation (14) containing a pair of complex-conjugated roots 
with a positive real part. At the same time, oscillatory type 
limit cycle SI such that phase difference cp varies within a 
limited range not exceeding 2n appears in phase space U. 
Cycle SI corresponds to a quasi-synchronous mode in the 
combined PLL and AGCL where periodic oscillations of phase 
variables are observed around equilibrium state Ai that has 
become unstable. Fig. 6a shows (<p,y)-projection of phase 
portrait and time waveform y(x) corresponding to the mode of 
cycle S 1 of model (3) for G(x)=5~Px. 



-2 66 




y 



-2.66 

-2.14 (p 3.14 O 

0.42 rn 1 0.42 

y 

-1 .72 

-314 (p 314 O 
1 09 I -it, I 1 OS 






-2 .32 

-3.14 (p 3.14 



y 

-2.82 




Fig. 6 Nonsynchronous modes of model (3) for G(x)=S-Px 
and y=0.86, u=0.727 (a); y=-0.54, u=1.6 (b), 2.63 (c) 
(ct=2, 0=1.1, 5=1.25, ei=l,s 2 =2) 
If, as a result of increasing u. or decreasing y, the system 
crosses line yp (where the saddle value ap=X4+Rexm<0, 
m= 1,2,3) a stable rotatory type limit cycle LI appears in the 
phase space U. In the mode of limit cycle LI phase difference 
cp rotates and variables y,z, and x periodically oscillate about 
certain mean values. 

Upon changing of parameters y and u. within the domain 
C s i (see Figs. 4, 5) limit cycle LI may disappear in the result of 
saddle-node bifurcation. After this bifurcation, the system 
undergoes a transition to the mode of rotatory chaotic attractor 
Wi via intermittence. Figs.6b, 6c show examples of (cp,y)- 
projections of phase portraits and dependences y(x) 
corresponding to the modes of attractors L, and W, of Model 
(3) for G(x)=5~Px. Modeling of System (3) indicates that 
complex multi turn-over rotatory type limit cycles exist in the 
phase space U when the system's parameters belong to domain 
C s i. Therefore, in domain C s i, combined PLL and AGCL have 
concurrent synchronous mode and asynchronous mode 
determined by stable equilibrium state and periodic or chaotic 
attractors of rotatory type. The initial conditions determine 
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which of these modes is realized in the system for parameters 
from the domain C s i. 

The quasi-synchronous mode of limit cycle SI exists for 
the values of parameters belonging to domain Cp restricted by 
the curves ys and yd (see Figs. 4, 5). Curve yd corresponds to 
the loss of stability of limit cycle S 1 as the result of the period- 
doubling bifurcation. The domain Cp is divided by bifurcation 
curve yp into domain Cp^ {max(ys,yp)<y<yd} and Cp 2 : 
{ys<y<min(yp, yd)}. For the parameter values from domain 
Cpi, limit cycle SI is globally stable. While for the parameters 
from domain Cp 2 -the Model (3) exhibits bistable behavior, the 
limit cycle SI and periodic or chaotic modes of attractors of 
rotational type concurrently exist in the phase space U. 
Transitions to the chaotic mode are realized via period- 
doubling bifurcations of rotatory limit cycles and so via 
intermittence. 

In domain C c , located to the right of line yd, the Model (3) 
exhibits complex dynamical behavior. For parameters from the 
domain C c quasi-synchronous and asynchronous modes of 
various complexities are realized in the combined PLL and 
AGCL. 

Let us track the evolution of the mode of limit cycle SI 
when parameter |i is increased. For this purpose, we use the 
results of numerical simulation of Model (3) for G(x)=o/x-[3x. 
The non-synchronous modes formed when parameter |i varies 
are illustrated by one -parameter bifurcation {|J,y} diagram of 
the Poincare mapping corresponding to rj=0.5, (3=1.1, a=2, 
8i=l, £2=2, y=0.1, and |ie(1.4;2.82) (Fig. 7), (cp,y) projections 
of the attractors' phase portraits, and fragments of dependences 

y(x)(Fig.8). 



1 .0 




-3.2 



2.07 ^ 2.8 

Fig. 7 Diagram showing evolution of quasi-synchronous mode of limit cycle 
S, during variation of \i for G(x)=a/x-px 




-3.14 (p 3.14 -314 (p 314 z 1000 



Fig. 8 Projections of phase portraits and fragments of time-realizations of 
oscillations corresponding to attractors of model (3) for G(x)=a/x~Px, u,=1.4 
(a); 1.95 (b,c); 2.0(d); 

2.07 (e,f); 2.21 (g,h); 2.55 (i); 2.62 (j); 2.82 (k,l) 

(a=0.5, (3=1.1, a=2, Ei=l, £s=2, y=0.1) 

The quasi-synchronous mode of limit cycle SI is the 
system's initial state. When |i increases, cycle SI exhibits 
period-doubling bifurcations, which finally end in a transition 
to the mode of oscillatory chaotic attractor V, (Figs. 8b, c). 
After that, the oscillations in the mode of attractor V! 
transform rigidly into the oscillations in the mode of 
oscillatory limit cycle Qj (Fig. 8d). When |i>2.07, the system 
switches from the mode of limit cycle Qj to the mode of 
chaotic attractor V 2 of oscillatory-rotatory type (Figs. 8e,f). 
Then, alternating chaotic and periodic asynchronous modes are 
observed; in the majority of the studied |i range, chaotic modes 
are realized. Examples of phase portraits and time realization 
y(x) of chaotic and periodic modes are shown in Figs. 8g-81. 

Numerical simulation of Model (3) reveals that this model 
may demonstrate such interesting dynamical phenomena as 
formation of attracting two-dimensional tori of oscillatory and 
rotatory types in the phase space U. These tori appear as the 
result of loss of oscillatory or rotatory limit cycles stability 
when a pair of complex-conjugated cycles, multiplicators 
crosses a unit circle. The tori correspond to two-frequency 
non-synchronous modes of combined PLL and AGCL. 

Fig. 9 illustrates the examples of the (cp,x)-projection of 
phase portraits, dependences x(x), and (y, x)-projections of 
Poincare mapping T generated by the phase trajectories of 
Model (3) corresponding to G(x)=5-Px and |3=1.1, 5=1.25, 
y=0.1, |i=5, £i=1.97, £ 2 =150. The results represented in Fig. 9 
characterize non-synchronous modes of the system formed 
when coupling parameter a varies. The quasi-synchronous 
mode of limit cycle S 2 (see Fig. 9a) is the system's starting 
state for oc=1.6. 

When a increases, the mode of oscillatory torus Tl (Fig. 
9b) appears from limit cycle S 2 . The phase portrait of mapping 
T (Fig. 9c) is characterized by the presence of stable closed 
invariant curve Tj. When 1.623<|i<1.92, the alternation of the 
mode of torus Tl and the modes of multi turn-over limit cycles 
of oscillatory type are observed. An example of nine-turn limit 
cycle S 3 is represented in Fig. 9d. Beginning with a=1.92, the 
distortion of the curve Ti is observed. This phenomenon 
indicates gradual transformation the mode of torus Tl to the 
mode of oscillatory type chaotic attractor P, (Fig. 9e). 
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Fig. 9 Attractors of model (3) for G(x)=5-|3x and 
a=1.6 (a); 1.65 (b,c); 1.9 (d); 1.946 (e); 2.08 (f,g,h); 
2.097 (i,j); 2.35 (k); 2.436 (l,m) 
(y=0.1, u=5, P=l.l, 5=1.25, el=1.97, e2=150) 
As a further increases, rigid transition of the system from 
the mode of chaotic attractor Pi to the mode of oscillatory 
eight-turn limit cycle S4 take place. With further increase of a 
the mode of cycle S 4 is transformed to the mode of chaotic 
attractor P 2 (Fig. 9f,g,h). Upon reaching the value a=2.088, 
rotatory phase trajectory scrolls in the attractor' s P 2 structure 
are appeared. This phenomenon indicates transformation the 
mode of oscillatory attractor P 2 to the mode of oscillatory- 
rotatory chaotic attractor W 2 (Fig. 9i, j). A still increase of a 
leads to the alternation of the modes of complex limit cycles of 
oscillatory-rotatory types (Fig. 9k) and the modes of chaotic 
attractors (Fig. 91, m) and following returning to the mode of 
periodic quasi-synchronous mode. 

The existence of rotatory tori in phase space U is revealed 
in the model (3) for G(x)=5~Px and P=l.l, 5=1.25, a=5, |i=5, 
£i=l, s 2 =150 when parameter y is varied. Figure 10 displaying 
the examples of (<p,y)-projection of phase portraits, 
dependences y(x), and (y,x)-projections of Poincare mapping 
Tcp of the plane cp=(p into the plane cp=cp +2jr produced by the 
trajectories of model (3) shows how mismatch y affects the 
process of the asynchronous mode of limit cycle L 2 (Fig. 10a) 
transformation during a decrease in y. 
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Fig. 10 Attractors of model (3) for G(x)=5-Px and 
7=0.853 (a); 0.845 (b); 0.819 (c,d); 0.818 (e,f,g); 
0.362 (h); 0.323 (i); 0.32 (j,k) 
(a=5, u=5, p=l.l, 5=1.25, £,=1.0, e 2 =150) 
At first, the mode of limit cycle L 2 is transformed to the 
mode of rotatory torus T2 (Fig. 10b). Torus T2 corresponds to 



asynchronous two-frequency mode of the combined PLL and 
AGCL. Further, the following phenomena are observed in the 
system with decreasing y: alternation of the mode of torus T2 
and the modes of asynchronous modes of multi turn-over 
rotatory limit cycles; gradual transformation of torus T2 to 
chaotic attractor W3 (Fig. 10c, d); alternation of chaotic (Fig. 
lOe, f, g) and periodic modes, the latter of which is determined 
by two-, four-, and three-turn rotatory limit cycles; soft 
transformation of the mode of three-turn rotatory limit cycle to 
the mode of rotatory torus T3; distortion of the torus T3 and 
following transition to the chaotic mode. A still decrease of y 
leads to transition of the system to the asynchronous mode of 
rotatory limit cycle L 3 (Fig. lOh) and following soft 
transformation of mode of cycle L 3 to the mode of rotatory 
torus T4 (Fig. lOi). When y continues to decrease, the mode of 
torus T4 is transformed to the mode of rotatory chaotic 
attractor W 4 (Fig. lOj, k). 

V. CONCLUSION 

In this paper, we have investigated the dynamical modes, 
bifurcation, and transitions to the chaotic behavior in a 
combined PLL and AGCL with the first- and second-order 
filters in PLL and the first-order filter in AGCL. Using the 
dynamical models of considered systems, we found that the 
systems exhibit a rich variety of dynamical modes including 
synchronous mode, periodic quasi-synchronous mode caused 
by the loss of the synchronous mode stability, asynchronous 
modes of rotatory and oscillatory-rotatory limit cycles, quasi- 
periodic quasi-synchronous and asynchronous modes 
corresponding to oscillatory and rotatory 2D tori in the phase 
space, chaotic modes that formed via period-doubling 
bifurcations, via intermittence, as well as rigidly formed 
through saddle-node bifurcations of limit cycles, via 
destruction of tori. These phenomena, in our opinion, are of 
importance for both basic and applied research of nonlinear 
dynamics of complex systems with phase control. It allows 
describing and explaining the behavior of the system when the 
synchronous state is cut off as a result of the system 
parameters perturbation. The wide variety of chaotic modes 
offers considerable possibilities of forming various frequency- 
and amplitude-modulated signals at the output of the system. 
Control the characteristics of generated signals can be easily 
realized by means of system's parameters. 
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